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Abstract 

Increasing evidence of a role of chronic inflammation in type 2 diabetes progression lias led to the development of 
therapies targeting the immune system. We develop a model of interleukin-1 (5 dynamics in order to explain principles of 
disease onset. The parameters in the model are derived from in vitro experiments and patient data. In the framework of this 
model, an IL-1 p switch is sufficient and necessary to account for type 2 diabetes onset. The model suggests that treatments 
targeting glucose bear the potential of stopping progression from pre-diabetes to overt type 2 diabetes. However, once in 
overt type 2 diabetes, these treatments have to be complemented by adjuvant anti-inflammatory therapies in order to stop 
or decelerate disease progression. Moreover, the model suggests that while glucose-lowering therapy needs to be 
continued all the way, dose and duration of the anti-inflammatory therapy needs to be specifically controlled. The model 
proposes a framework for the discussion of clinical trial outcomes. 
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Introduction 

Despite more than 350 million patients worldwide and the 
concomitant expensive socioeconomic burden, the pathogenesis of 
type 2 diabetes (T2D) is not yet completely understood. T2D is a 
progressive disease. The most important physiological components 
of T2D are insulin resistance, which is characterized by impaired 
response to insulin in insulin-sensitive tissues, and P-cell failure, 
which is characterized by P-ceU dysfunction and reduced P-ceU 
mass. The progression of T2D is clearly divided into at least two 
phases, pre-diabetes and overt diabetes [1-5]. In the pre-diabetes 
phase, insulin resistance is compensated by increased single P-ceU 
secretion capacity and/ or P-ceU number. If insulin resistance is not 
completely compensated, the blood glucose level would grow 
slowly, manifested as higher fasting glucose (impaired fasting 
glucose, IFG) and/or higher post-load glucose (impaired glucose 
tolerance, IGT) [6]. Overt T2D is characterized by compensation 
failure and continuous loss of functional P-ceUs [7,8], hence 
accompanied by continuously aggravated hyperglycaemia. 

Although insulin resistance is usually present in the early phase 
of pre-diabetes, it is the pace of P-ceU faUure that determines the 
onset of overt T2D [9] . The mechanisms leading to the transition 
from pre-diabetes to overt T2D are unclear [10]. However, there 
is evidence that the transition from P-ceU compensation to P-ceU 
faUure happens in a comparably short time span [2,4], typicaUy 
within 3 years [2]. This is further supported by a recent 



longitudinal study in a large population [5]. The trajectory of 
glycaemia before diagnosis of T2D was shown to be composed of a 
slow and stable adaptation, which lasts 12 years, foUowed by a 
rapid rise of glucose to overt T2D within 2 years [5] . Once overt 
T2D is started, hyperglycaemia continues to worsen regardless of 
treatments based on oral anti-diabetic agents [11-13]. 

The evidence of a role of inflammatory responses in the 
pathogenesis of T2D was increasing in recent years. Interleukin- 1 P 
(IL-ip) has been reported to contribute to P-ceU faUure [14—18]. P- 
ceUs themselves secrete IL-ip upon glucose stimulation [14]. 
Furthermore, IL-ip stimulates its own production in P-ceUs [17] 
and attracts macrophages [18] which can act as an extra source of 
IL-ip and other cytokines. Although it is currently unclear 
whether inflammatory responses are a primary cause or a 
secondary effect in T2D progression, therapies targeting IL-ip 
have shown encouraging progress albeit diverse results in different 
clinical trials [19-26]. 

These results motivated our working hypothesis that the pre- 
diabetic and overt T2D might be characterised by two quaUta- 
tively different states and that IL-ip is a potential candidate for 
promoting the transition between these two states. The hypothesis 
is consistent with results from clinic trials in which IL- 1 P blockade 
by interleukin- 1 -receptor antagonist (IL-lRa), a naturaUy occur- 
ring competitive inhibitor of IL- 1 P, induced sustained improve- 
ments of P-ceU function and the systemic inflammation state in 
patients with a mean disease duration of 1 1 years, even 39 weeks 
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Author Summary 

Insulin resistance and relative insulin deficiency are two 
hallmarks of type 2 diabetes. While insulin resistance is 
always present in the early phase, it is p-cell failure that 
determines the pace of the disease onset. Increasing 
evidence that the immune system is activated and plays an 
important role in type 2 diabetes has stimulated efforts of 
developing drugs targeting inflammatory cytokines. We 
built a model to describe the principles of type 2 diabetes 
onset under the influence of interleukin-1 p. The disease 
onset is understood in terms of bifurcation. It is found that 
inflammatory cytokines are required to be suppressed for a 
limited time only, while glucose has to be controlled over 
the long term. These structural insights may serve as a 
guideline for future clinical trials. 

after cessation of treatment [19,20]. IL-lRa competes with IL-ip 
for IL- 1 receptors but does not trigger any signalling event. It is an 
important regulator of the effect of IL- 1 P in many cell types, 
including human pancreatic P-ceUs which secrete IL-lRa them- 
selves [27,28]. In several liver, autoimmune, and infectious 
diseases IL-lRa is a better indicator of disease severity than IL- 
IP [29]. Results from a recent longitudinal study show that an 
accelerated increase in circulating IL-lRa starts about 5 years 
before diagnosis of T2D [30], coinciding with the accelerated 
deterioration of insulin-sensitivity and the compensation of P-ceU 
function [5] . In addition, long-term effects of temporary intensive 
insulin therapy on diabetes remission (or at least a prolonged 
normoglycaemia phase) in newly diagnosed diabetic patients may 
be partially attributed to the anti-inflammatory activity of insulin 
[3 1] . AH these results highlight the role of inflammatory responses 
in the pre-diabetes phase and the possibility that the long-term 
effects induced by temporary anti-inflammatory therapy are 
mediated by a switch, i.e. a sudden transition, between the 
qualitatively different compensation and overt T2D states. 
Interestingly, in vitro experiments show that the effect of IL- 1 P 
on P-ceU mass and insulin secretion is twofold [32]: low 
concentration of IL- 1 P stimulates P-ceU proliferation, inhibits P- 
ceU apoptosis and enhances glucose-stimulated insulin secretion 
while high concentration of IL-ip has opposite effects. In other 
words, IL-ip may contribute to both the P-ceU compensation 
phase and the P-ceU failure phase. 

Here, we present a mathematical model showing the 
possibility that T2D onset is induced by a sudden transition of 
IL-ip to a high level. The IL-ip switch results from the 
coexistence of two different states, which is often termed as bi- 
stability. Bi-stable switch models have been widely used in 
modeling developmental processes, such as cell cycle progression, 
cellular differentiation and apoptosis. The characteristics of bi- 
stable switches include sudden transition, and hysteresis. The 
latter means that the state of a system depends not only on the 
current environment but also on its past state. We combine the 
IL-ip bi-stable switch model with a previously published T2D 
progression model [33], by assuming extrapolated P-ceU turnover 
rates caused by exogenous IL-ip [32]. Similar concepts 
describing diabetes progression by different steady states have 
been already developed for both type 1 and 2 diabetes [33-35]. 
Our model falls into the same class as those previously published 
by Topp et al. [33] and De Gaetano et al. [35], which considered 
the evolution of P-ceU mass, insulinemia and glycaemia over a 
time-scale of years and rely on glucose toxicity. P-ceU mass is 
controlled by an empirical parabolic function of glucose in [33]. 
Subsequently a physiologicaUy amenable concept of pancreatic 



reserve was introduced, which controls the direction of P-ceU 
mass change [35]. Both P-ceU mass and the pancreatic reserve 
were modeled phenomenologically as functions of glucose, where 
some parameters bear the potential to act as bifurcation 
parameters. The model presented here is buUt on these results. 
The phenomenological implementation of glucose toxicity is 
replaced by a model of IL- 1 P, which is explicitly defined and 
relatively easy to measure. We focus on an IL-1 P bi-stable switch 
giving rise to hysteresis which is not discussed in previous studies. 
IL-ip hysteresis turns out to be of utmost importance because it 
widely shapes the strategy of anti-inflammatory therapies. While 
the complexity of the disease wiU never be captured by a 
mathematical model, the structural insights suggested by the 
model will help exploiting treatment approaches more efficiently 
and elaborating the effect of different therapies on P-ceU mass 
and disease progression. 

The IL- 1 P bi-stable switch model implies that T2D onset and 
later irreversible progression of P-ceU failure arise from altered 
stability properties of the P-ceU compensation state. Therefore, 
the strategy of glucose control would be more effective in the 
compensation phase than after the onset of overt T2D. This 
model result is consistent with recent data from the Diabetes 
Prevention Program Outcome Study (DPPOS) showing that 
IGT people who once regressed to normoglycaemia have a 56% 
lower relative risk of developing diabetes [36]. While glucose- 
lowering therapy bears the potential of stopping disease 
progression in the compensation phase, the IL-ip bi-stable 
switch model suggests that a combined glucose-lowering and 
anti-inflammatory therapy is necessary to stop the loss of 
functional P-ceU mass as well as hyperglycaemia progression in 
the overt T2D phase, and eventually reestablish the compen- 
sation phase. 

Materials and Methods 

A full description of the model (Fig. 1), the assumptions, and the 
derivation of aU model parameters from experimental data 
(summarized in Table 1) are provided in this section. Scripts of 
aU calculations are provided in the supplement, including the 
model fUe and the parameter estimation procedures. 

Overview 

The mathematical model (Fig. 1) is constructed in three steps. 
First, we build the core of the model which describes the 
competition between IL- 1 P and IL- 1 Ra, where glucose stimulates 
IL-ip production. Since the effect of IL-1 P on P-ceUs is mediated 
by the binding to its receptor, which is also the target of IL-lRa 
antagonism, the fraction of IL- 1 P bound receptor F is used as a 
measure for IL-ip stimulation. 

Second, the interaction between glucose (G) and insulin (I) is 
modeled with P-ceU mass being a parameter. This subsystem is 
mainly adopted from [33], which is a single-compartment 
model justified by the slow dynamics of glucose and insulin over 
a time-scale of days to years. G and I are associated with long- 
term fasting glucose and insulin levels. As short-term (daily) 
fluctuations are averaged out, the model cannot describe 
disease onset derived from details of glucose and insulin 
dynamics. 

Third, the two subsystems which describe the inflammatory 
signals and glucose-dependent insulin secretion are connected. For 
that purpose the IL- 1 P-dependent P-ceU turnover rates are 
adopted from human islets which were cultured with exogenous 
IL- 1 p. The data suggest a bimodal effect of IL- 1 P on the P-ceU 
mass [32]: IL-ip stimulates P-ceU proliferation and inhibits P-ceU 
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Figure 1. Model scheme: IL-IRa {A) and IL-ip (Z.) compete for IL-1 receptors giving rise to a fraction of IL-ip bound receptors [f) 
which determines subsequent signalling in p-cells (B). These control insulin (/) release and, via the influence of insulin resistance, blood glucose 
level (G). Arrows: activation effect, line with bar end: inhibition effect. Kinetic terms corresponding to each interaction are labelled. Variables evolving 
on different time scales are marked by different colours. 
doi:1 0.1 371 /journal.pcbi.1 003798.g001 



apoptosis when presented in low concentrations. Conversely, at 
high concentrations IL-ip enhances P-cell apoptosis and reduces 
P-cell proliferation. The unknown model parameters are deter- 
mined from steady state conditions of the IL- 1 p/IL- 1 Ra 
subsystem and from the natural history of fasting glucose in the 
pre-diabetic phase [5]. 

There are three different time scales in the model (Fig. 1). The 
IL- 1 p/IL- 1 Ra subsystem evolves on the fastest scale, since ligands 
and receptors bind/unbind on a time scale of seconds and proteins 
are synthesised/secreted on a time scale half an hour. The 
glucose/insulin subsystem evolves on an intermediate time scale of 
days [33]. P-ceU evolves on the slowest scale. The turnover rate of 
P-ceU is rather slow: typically, one or two cells divide per month. 
Since our main interest lies on long-term evolution of fasting 
glucose, the IL-ip/IL-lRa subsystem is assumed in steady-state 
equilibrium and P-ceU mass is assumed to be constant on this 
relatively fast time scale of fasting glucose. The different time- 



scales in the model make it possible to separate the model into 
different subsystems and justify our three-step modelling approach. 

The antagonist effect of IL-IRa {A) 

All effects of IL-ip (L) are mediated by the IL-1 receptor 
and are blocked by IL-IRa (A). IL-IRa competes with IL-ip 
for the receptor but does not trigger any signalling event. 
Therefore, ligand-receptor binding and unbinding are consid- 
ered to describe the antagonist effect of IL-IRa. The ligand- 
receptor binding/unbinding happens on a time scale of 
seconds and is fast compared to other time-scales in the 
model, such as IL-ip and IL-IRa production and secretion 
which happens on a time-scale of tens of minutes. Therefore, it 
is reasonable to assume that ligand and receptor are always in 
equilibrium, which enables us to describe the fraction of IL-1 P 
bound receptors (F) by 
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Table 1. Summary of parameters used in the model. 



Description 


Symbol 


Value 


Unit 


Ref. 


glucose 


6 








initial value 




5.47(A) 


mM 


[5] 






5.25(W) 


mM 




glucose production 


^0 


48 


mM d"' 


[33] 


insulin-independent Glucose uptake 




1.44 


d-' 


[33] 


insulin / 


initial value 




70.73(ft) 


pM 


DBSS 






74.28(W) 


pM 




maximal rate of insulin secretion by p-cells 


h 


300.02 


pM d"' 


[33] 


insulin clearance rate 


d, 


432 


d-^ 


[33] 


glucose stimulated insulin secretion 


s(G) 










Kg 


7.86 


mM^ 


[33] 


p-cell mass B 


initial value 




312.1(ft) 


mgmg 


DBSS 






346.7(W) 




DBSS 


ratio of proliferation to apoptosis 


X 


2.33 




DBSS 


p-cell proliferation 


P(F) 










Op 


-0.048 




by fit 






-0.46 




by fit 






0.62 




by fit 


p-cell apoptosis rate 














0.047 




by fit 




ba 


0.57 




by fit 




Ca 


2.36 




by fit 


p-cell turnover time constant 


T 


0.00052 


d-' 


FBNH 


IL-IRa 


A 








initial value 




6000(A) 


pg mP^ 








5982(W) 


pg mP^ 


DBSS 


IL-lp stimulated IL-IRa production 




k, 


9.47e5 


pg mr' d"' 


DBSS 






3.1 2e7 


pg mr^ d^' 


DBSS 






2.65e4 


d-' 


DBSS 


degradation rate of IL-IRa 




166.36 


d-' 


[29] 


dissociation constant of IL-IRa 




1300 


pg mP' 


[29] 


IL-ip 


L 








initial value 




12(A) 


pg mP' 


* 






11.25(W-) 


pg mP^ 


DBSS 


inhibitory effects of basal IL-1 p mRNA 


b{L) 










u 


5.8 




by fit 




r 


2 




by fit 


coefficient relating IL-1 p mRNA level to its concentration 


k. 


0.047 


pg"^ ml 


DBSS 


glucose stimulated IL-lp production 


g(G) 










s 


2.77 




by fit 




V 


3.52 


mM 


by fit 




ks 


3.745e3 


pgmP' d"' 


DBSS 


IL-lp auto-stimulation 


m 










f 


3 




by fit 




Kf 


0.1575 




by fit 




ke 


1 .048e9 


pgmP' d"' 


DBSS 


degradation rate of IL-lp 


dt 


55.45 


d-' 


[63] 
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Table 1. Cont. 





Description 


Symbol 


Value 


Unit 


Ref. 


dissociation constant of IL-1 p 


Kl 


1300 


pg mT^ 


[291 


insulin resistance 


R 










Ro 


9.643 


pM d 


[33] 




m 


1.688 


pM d 


FBNH 



Parameter values of the mathematical model and references. A: initial value of incident diabetes cases; W: initial value of non-diabetic control; FBNH: fitted by glucose 

history; DBSS: determined by steady state conditions; 

*: parameter discussed in the section "steady state conditions". 

doi:l 0.1 371/journal.pcbi.1 003798.t001 



(1) 



where Ki^ and Kp^ are the dissociation constants of IL-1 (5 and 
IL-lRa, respectively. We assume a constant number of IL-1 
receptors on the P-cell membrane. Consequently, F can be 
considered as a measure of IL-1 (5 stimulation. 

Estimate of endogenous IL-ip [L) and IL-IRa [A] in 
cultured islets 

F is used in the following to fit certain functions in the model to 
corresponding experimental data. In these experiments, dififerent 
amounts of exogenous IL- 1 P were added to cultured human islets 
and the corresponding effects were measured, while endogenous 
IL-IRa and IL-ip were not measured. It is safe to neglect 
endogenous IL-ip and IL-IRa when exogenous IL-ip is 
dominant, which is the case for most measured data points. 
However, it is necessary to estimate the endogenous IL- 1 P and IL- 
lRa concentrations for data points with zero exogenous IL-ip. 

Secretion of IL- 1 Ra from cultured human islets was measured 
after 4 days which led to the endogenous IL- 1 Ra concentration of 
114 pg/ml [32]. Endogenous IL-ip is estimated to be 0.228 pg/ 
ml using the ratio of IL-lRa/IL-1 P of 500 [27]. 



stimulation [17] (&(L)). IL-ip degrades with a rate di^. In the 
following, the functions g[G), /(F), and 6(L) are derived from 
experimental data. The parameters and are determined in 
the Section "steady state conditions". 

1. IL-lp induction by glucose (g{G)). We use a Hill function 
g{G) 



g{G) = ag 



(4) 



to fit data from cultured human islet describing the increase of IL- 1 P 
mRNA (Fig. 4A in [17]) as a function of glucose (Fig. 2 A). This 
implicidy assumes that IL- 1 P mRNA level is proportional to IL- 1 P 
concentration. Since the available data do not include glucose levels 
below 5.5 mM, we assume that half of the control rate is achieved at 
3 mM (marked in Fig. 2A by the black square), which is the usual 
limit of hypoglycaemia. The choice of this value does not alter the 
model results. Other values, such as 2 or 4 mM, only affect the 
values of the model parameters (Ai-e) determined by steady state 
conditions (see Section below). We ignore the data point at 40 mM 
glucose, since such a high glucose level is not physiological, = 1 is 
used in the model, since the constant is merged with kf. 

2. IL-lp auto-stimulation [1{F)). IL-ip auto-stimulation is 
described by 



IL-IRa [A) dynamics 

IL-lRa {A) stimulation by IL-ip (L) is bimodal [32]. IL-IRa is 
induced by IL-ip at exogenous concentrations below 20 pg/ml. 
At higher concentrations, IL-ip restores IL-IRa to nearly the 
basal level. This nonlinear effect is phenomenologically captured 
by the nonlinear equation 



AA 
dt 



--ki^kiF-k^F^A- cIaA 



(2) 



where (/a is the natural degradation rate of IL-IRa. The 
parameters ki to A3 are determined in the Section "steady states 
conditions" (see below). 

IL-1 p (/.) dynamics 

IL- 1 P dynamics are described by 



dL 



--biL)(ksg{G) + kel{F))-di^L 



(3) 



Glucose induces P-cells to secrete IL- 1 P [17] (A5 g{G)). IL- 1 P is also 
secreted from islets via auto-stimulation [1 7] (kf, 1(F)). High levels of 
IL- 1 P mRNA inhibit the stimulating effect of glucose and of auto- 



1(F) = a, 



Kf' 



(5) 



and is fitted to data (Fig. 4B in [17]) which determine the increase 
of IL-ip mRNA of cultured human islets in dependence on 
exogenous IL-ip concentrations (Fig. 2 B). a;= 1 is used in the 
model, since the constant is merged with kf,. 

3. Inhibition of the effects of glucose and IL-ip by IL-lp 
mRNA [b[L)). The stimulatory effect of glucose is inhibited by basal 
IL- 1 P mRNA [17]. IL- 1 P auto-stimulation Ls linearly related to glucose 
stimulated IL-ip in human islet preparations [17]. This linear 
relationship suggests a common regulatory mechanism of glucose- 
stimulated and auto-stimulated IL-ip production. This is represented 
in Eq. (3) by basal IL-1 P niRNA-dependent inhibition of both glucose 
stimulation and IL- 1 P auto-stimulation using the sigmoidal fimction 



b(L)-- 



at 



u + (k4LY 



(6) 



where the scaling parameter A4 relates IL- 1 P mRNA level to IL- 1 P 
protein concentration. Eq. (6) is fitted to the data (Fig. IE in [17]) as a 
function of basal IL- 1 P (Fig. 2 C) and A4 is determined in the Section 
"steady state conditions", fit/, = 1 is used in the model, since the constant 
is merged with kg and kf,. 



PLOS Computational Biology | www.ploscompbiol.org 



5 



August 2014 | Volume 10 | Issue 8 | el 003798 



Modelling Pre-diabetes Progression 



A 
2 

1.5 



o 

1 0.5 



E 
2.5 



proposed data 
data 

fitted curve 



10 20 30 

glucose (mM) 



40 




2 4 6 

basal IL-ip mRNA 




1.5 



0.5' ■ ■ ■ ■ 1 

-8 -6 -4 -2 0 

log(F) 

Figure 2. Parameter fit from experimental data. A: Data are reproduced from Fig. 4A in [17]. The function g(G) in Eq. (4) is fitted to these data. 
Note a = 1 is used in the model, since the constant could be merged with k^. B: Data were reproduced from Fig. 4B in [1 7]. The function /(F) in Eq. (5) is 
fitted to these data. Note a = 1 is used in the model. C: Data are reproduced from Fig. 1 E in [1 7]. Eq. (6) is fitted to these data. Note the interchanged 
axes and a = 1 is used in the model. D: Data are reproduced from Fig. 1C in [32]. The function a(F) in Eq. (13) is fitted to these data. E: Data are 
reproduced from Fig. IB in [32]. The function p(F) in Eq. (14) is fitted to these data. 
doi:1 0.1 371 /journal.pcbi.1 003798.g002 



Glucose (G) and insulin (/) dynamics 

G and I are quantities associated witli long-term changes (days to 
years) of fasting levels of glucose and insulin, respectively. Therefore, 
a single-compartment model for glucose and insulin is used. The 
ordinary differential equations describing glucose and insulin are 



adapted from Topp et al. [33] with the modification that the 
dynamics of insulin resistance {R) is included as input to the model. 

1. Glucose dynamics (G). Glucose is produced with a 
constant rate (ro) and consumed by an insulin-dependent {I/R) and 
an insulin-independent rate {do) 
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dG 

lit 



--ro- 



(7) 



where R denotes insulin resistance. 

2. Insulin dynamics (7). Insulin dynamics is modeled by a 
mass balance equation 



d/ 
dt ^ 



--hs{G)B-diI 



(8) 



with a secretion rate s(G) and a clearance rate (di). Secretion is 
stimulated by glucose according to the sigmoidal function 



.5(G) = y- 



G2 



(9) 



and is proportional to the P-cell mass (5). The maximal rate of 
insulin secretion h, the insulin clearance rate di, and the half 
concentration Kq of glucose stimulated insulin secretion are taken 
from [33]. 

3. Insulin resistance {R). Insulin sensitivity provided in [5] 
is homeostatic model assessment (HOMA) index. The one 
compartment model of glucose and insulin dynamics is a 
simplification of the minimal model [37]. Insulin resistance 
derived from the minimal model exhibits a linear relationship to 
the HOMA-IR index on logarithmic scale [38]. Hence, 



R = Ro 



Sh.o 



(11) 



is used to convert the HOMA index value into R. Here S// is the 
insulin sensitivity index (HOMA2 %S) reported in [5]. The 
parameter m is determined by fitting the model to the natural 
history of fasting glucose in the pre-diabetic phase [5] . 

jS-cell mass dynamics (6) 

The dynamics of the P-cell mass 



dB 
dt 



T%(F)-/fl(F)] 



(12) 



is controlled by IL-lfi induced fi-ceU apoptosis and proliferation. 
The apoptosis rate a(Fj 



a(F) = aJog'{F) + b,log(F) + c, 



(13) 



and the proliferation rate p{F) 



p(F) = flp log2(F) + bp log(F) + Cp 



(14) 



are fitted to data from human islets [32], which were cultured for 4 
days with different IL- 1 P concentrations applied at day 0 (Fig. 2 D 
and E, respectively). 

We assume that the initiation of apoptosis and proliferation is 
faster than IL- 1 P-induced IL-ip and IL-lRa secretion. Thus, the 
IL-ip and IL-lRa dynamics are de-coupled from P-cell turnover 
and the initial F (at day 0 of the experiment) can be used as an 
approximation in fitting the data. Since the apoptosis and the 
proliferation rates depend on log(F), eventual modifications of F 
caused by newly secreted molecules would anyway remain small 
on the logarithmic scale. 

P-cell apoptosis and proliferation were measured as percentage 
of control (i.e. no IL-ip added) [32] such that the absolute rates 
remained undetermined. Scaling of the rates is represented by an 
extra parameter T (Eq. ( 1 2)), which is fitted to the natural history of 
fasting glucose in the pre-diabetic phase [5]. The relative 
importance of apoptosis and proliferation is controlled by the 
parameter X in Eq. (12), which is determined in the Section 
"steady state conditions". 

Steady state conditions 

While a large number of model parameters could be determined 
by experimental constraints, the parameters Ai_g and A remained 
undetermined and are derived from steady state conditions for 
Eqs. (2,3,12). Eq. (12) is used to determine /.. Eq. (2) and (3) each 
exhibit two stable steady states and a bifurcation point, such that 
these 6 additional conditions can be used to determine Ai_b. 

1. Steady state of the compensation phase. Reports from 
the Whitehall II longitudinal study [5] show that 1 3 years before 
diagnosis of T2D, the fasting glucose level of the incident diabetic 
participants is already higher than that of control (5.47 vs. 
5.25 mM, Table 2). The corresponding serum IL-lRa concentra- 
tion of the incident diabetic participants is 300 pg/ml [30]. Using 
the reported in vitro ratio of IL-lRa/IL-ip, which is 500 [27], the 
corresponding serum IL-ip concentration is estimated to be 
0.6 pg/ml. This value is consistent with the average serum IL-ip 
level (0.57±0.93 pg/ml) reported in [39]. Since IL-ip/IL-lRa are 
locally produced in the pancreas islets and their half-lives are short 
(5~10 minutes), we assume 20 fold more IL-ip/IL-lRa in the 
pancreas than in the blood under healthy conditions (note that the 
blood-tissue ratio might change after T2D onset), inferring 1 2 and 
6000 pg/ml (Table 2) of IL-ip and IL-lRa, respectively, in the 
islets. As this choice is arbitrary, other choices of the tissue-blood 
ratio, such as 200 or 2000, with fixed ratio of IL-lRa/IL-1 P, were 
tested. The ratio does not change the model behaviour since the 
fraction of IL- 1 P bound receptors F rather than the level of IL- 1 P 
determines the inflammatory stimulus. However, the choice of the 



Table 2. Steady state conditions for the IL-ip/IL-1 Ra 


subsystem. 








Glucose(mM) 


IL-ip (pg/ml) 


IL-1 Ra (pg/ml) 




(G) 


(L) 


(A) 


compensation 


5.47 


12 


6000 


bifurcation 


5.84 


Ltj (unknown) 


At, (unl<nown) 


overt diabetes 


13.4 


12*55.8 


6000*1.8 


doi:10.1371/jou 


rnal.pcbi.1003798.t002 
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tissue-blood ratio affects tlie values of the model parameters {X and 
derived from steady state conditions. 

2. Steady state of P-cell mass determines i. P-ceU 
proliferation and apoptosis are assumed to be exacdy balanced 
at a healthy glucose level (4.5 mM), which determines the 
parameter 1 assuming Eq. (12) in steady state using 
f o(Go - 4.5 mM). If the glucose level is less than 4.5 mM, the 
apoptosis rate would surpass the proliferation rate which would 
decrease the P-ceU mass. When the glucose level increases, more 
IL-ip is produced which stimulates P-ceU proliferation and 
increases the P-ceU mass. Increasing insulin resistance detunes 
this equilibrium. 

3. Steady state at the critical (bifurcation) point. A 
sudden rapid increase in fasting glucose happens about two 
years before diagnosis of T2D. The average fasting glucose level 
at the inflection point of glucose history is 5.84 niM [5] 
(Table 2). 

4. Steady state in the overt diabetes phase. IL-1 P mRNA 
expression was measured in P-cells from T2D and non-diabetic 
cadaver donors [17]. Real-time quantitative PGR revealed in 
average 55.8 fold more IL-ip mRNA in T2D patients compared 
to healthy controls. As discussed above, we assume a linear 
relation of mRNA and protein concentration. The average fasting 
glucose of those diabetic donors was 13.4 mM (Table 2) [17]. 
Because of the bimodal dependency of IL-lRa on IL-ip (see 
Fig. 2E in [32]), two IL-lRa levels may be consistent with the 55.8 
fold increase of IL-ip: Either the basal level of IL-lRa or a 1.8 fold 
increased level. We tested both solutions by parameter scanning 
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(see below) and only the latter solution yielded consistent results 
and is, therefore, taken in the model (Table 2). 

5. Steady state equations. The steady state equations for 

Eq. (2) and (3) read 

0 = ki +k2Fi — kjFi^A[—df^Ai 
0 = b(U){k5g(Gi) + k6l{F,)) - di^Li 

where i= {c,b,d} correspond to the compensation, bifurcation and 
diabetes state, respectively. The values are summarized in Table 2. 
Note that IL-ip (L],) and IL-lRa (A]-) are not determined by 
experimental data. 

6. Determine the range of IL-ip and IL-lRa at bifurcation 
by parameter scan. The parameter plane is scanned in 
order to identify possible IE- 1 P and IE- 1 Ra concentrations at the 
bifurcation point. The values are restricted by the side condition that 
a switch of IL-1 P and IL-lRa exists (Fig. 3 A and B). This is done by 
numerically checking the determinant of the Jacobian and then 
checking the number of solutions on both sides of the bifitrcation 
point Gf, (the Madab script is included in Dataset SI). While a switch 
of IL-1 P exists in a large area of the L^-A], plane (red region. Fig. 3 
A), the required switch to a 55.8 fold IL-1 P level is only achieved for 
a small subset of values (Fig. 3 B). Most parameter combinations led 
to a smaller switch of IE- 1 p. This determined the choice of the values 
for L|, and^i,. Note that the scan for dififerent basal lE-ip/IL-lRa 
levels ( 1 20 /6 X 1 0* or 1 200/ 6 X 1 with IE- 1 Ra/IL- 1 P ratio frxed at 
the reported value [27]) yielded similar results (Fig. 3). 




Figure 3. Parameter determined by scanning. A: Random points are selected in the /.b-zAb plane (normalised with and /AJ and used in the 
steady state equations (Eq. (2,3)) to determine k^_^. If k^_^ are non-negative, points either induce an IL-1 (3 switch (red) or do not (blue). B: The IL-1 13 
level /-d/i-c 3fter the switch is shown as a function of IL-lRa concentration at the bifurcation point. 
doi:1 0.1 371/journal.pcbi.1 003798.g003 
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Figure 4. Fitting the model to the fasting glucose history to determine m and r. Shown are the best fit results of incident diabetes cases 
before T2D diagnosis and the data (black line) (red line cross) and non-diabetic controls (red line) (A), interpolated insulin resistance as model input 
(B), corresponding behaviour of p-cell mass (C) and IL-1 p (D). 
doi:10.1371/journal.pcbi.1003798.g004 



Fitting the natural history of fasting glucose 

The last two undetermined parameters, m and T defined in Eq. 
11 and 12, are determined by fitting the model to the fasting 
glucose history of both incident diabetes cases and the non- 
diabetics controls [5]. The best fit is shown in Fig. 4A. 
Interpolated insulin sensitivity [5] (Fig. 4B) which is an input to 
the model during fitting, the behaviour of the P-cell mass (Fig. 4C) 
and of IL- 1 P (Fig. 4D) are also shown. The fitting is implemented 
in the SBTOOLBOX2 [40] using a differential evolution 
algorithm, which reports multiple possible parameter sets 
(Fig. 5A). A strong correlation between m and T was identified 
(Fig. 5B), which gives rise to large variations of glucose and P-ceU 
mass in the simulations of overt diabetes (Fig. SI). In other words, 
the pre-diabetes glucose trajectory does not allow for a precise 
determination of ?« and T. The values from the best fit are reported 
in Table 1 and used throughout the paper. 

A complementary approach 

The above approach demonstrates that the bifurcation of the IL- 
1 p/IL- 1 Ra subsystem is sufficient to account for the fitted data set. To 
investigate whether the bifurcation of the IL-1 p/IL-lRa subsystem is 



necessary, a complementary approach is employed: the model is fitted 
to the data without bifurcation as a presupposed condition. The 
steady state equations for Eq. (2) and (3), in the compensation phase 
and in the overt diabetes phase, allow reducing the degree of freedom 
from 8 {ki-(j, m and t) to 4: A;, A2, k;, and are expHcitiy expressed 
as functions of and k4, and only k;, k4, m and T are fitted. Strong 
correlations are found to exist between m and T, which cause large 
variations of glucose and P-cell mass in the simulations of the overt 
diabetes (Fig. S2 A and B). Strong correlations are also found between 
711 3iidk4, which cause large variations of the time point of the IL-ip/ 
IL-lRa switch (Fig. S2 C and D). Please note that the complementary 
approach found, although with low probability, cases allowing an 
extra transition (see Fig. S2 C and D). 

Although some parameters cannot be determined precisely, the 
key qualitative property, i.e. the switch of the IL-ip/IL-lRa 
subsystem before T2D diagnosis, remains robust. Taking together 
the results of both complementary approaches, it is suggested that 
the bifurcation of the IL-ip/IL-lRa subsystem is both sufficient 
and necessary to account for the data set, in the framework of the 
current model. Please note that despite the model being vahdated 
by quantitative experimental data, the results presented here 
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Figure 5. There exist multiple reasonable fits (shadow) of incident diabetes cases before T2D diagnosis and the data (black line) 

and non-diabetic controls (red line) (A). Fitted parameters exhibit correlations (B). 

doi:10.1371/journal.pcbi.1003798.g005 



remain of qualitative nature. In particular, in the overt T2D state 
we expect that additional physiological process as well as 
medication will alter the results on a quantitative level. 

Results 

The IL-ip/IL-1Ra subsystem has multiple steady states 

An IL-ip switch is observed in a large range of parameter 
values including the physiological relevant range. Below the 
critical glucose level of 5.84 mM the model gives rise to five 
steady states, three stable and two unstable ones (Fig. 6A). The 
stable steady state with lowest IT- 1 P (green in Fig. 6 A) is stable 
only if glucose is below the critical level. It is interpreted as 
corresponding to the compensation phase. The hyperglycaemia 
parts of the two other stable steady states may be associated with 
overt T2D. 

NuUcline analysis facilitates understanding the transition between 
the different inflammatory levels. Steady states appear as the 
crossing points of tlie IL-lRa versus IL-ip nuUclines where both 
quantities become static (Fig. 6B). The IL-ip nuUcHne (Fig. 6B red 
line) exhibits a nadir followed by a peak, which is the typical feature 
of an auto-stimulation species. The IL-lRa nuUchne (Fig. 6B black 
line) exhibits one peak. This implies that at large IL-ip its further 
increase is accompanied by decreased IL-lRa. The measured non- 
linear dependence of IL-lRa on IL-ip [32] is at the origin of the 
third stable steady state with highest IL- 1 p. If IL- 1 Ra is stimulated 
by IL-ip in a linear manner, only two stable steady states would 
have been found. Such a dysregulation of IL-lRa may also happen 
in other diseases involving IL-ip, such as chronic myelogenous 
leukaemia and hairy cell leukaemia [29]. Hence, dysregulation of 
IL- 1 Ra at high IL- 1 P concentrations might be associated with a 
boost of disease progression. 

The compensation steady state is lost upon progression 
to overt T2D 

The stable steady state with lowest IL- 1 P corresponds to the pre- 
diabetic state. As glucose influences IL- 1 P production [14,17] the nadir 
of the IL-ip nuUcHne moves up in response to increased glucose 



(Fig. 6B inset) while the IL-lRa nuUcMne remains unchanged (Fig. 6B, 
black line). The two involved steady states approach each other and 
ultimately merge into one steady state (Fig. 6B inset, blue Mne). At even 
higher glucose this steady state disappears and the system makes a 
transition (Fig. 6B, black arrow) to the next available stable steady state, 
which is characterised by substantially higher IL-ip and moderately 
increased IL-lRa. This transition is interpreted to correspond to the 
sudden transition from the compensation phase to disease progression. 

At even higher levels of IL-ip, the effect of glucose on IL-ip 
production is inhibited by IL-ip itself (represented in Eq. 3 by h(l)), 
such that the IL- 1 P nullcHne remains unaffected by glucose at high 
IL- 1 P concentrations. Therefore, the current model only reflects the 
transition from the compensation phase to the mild disease phase and 
not to the strong disease phase. An extended model including detailed 
IL- 1 P mRNA dynamics and other factors like free fatty acids, leptin 
or tissue resident macrophages might generate two transitions. 

Glucose-lowering therapy cannot cure or stop 
progression of T2D 

Glucose-lowering therapy, simulated in silico by fixing the level of 
glucose at normal level (5.47 mM), is represented by the IL-ip 
nuUcHne which now, as in the pre-diabetic state, intersects with the 
IL- 1 Ra nullcline (Fig. 7 A) around its nadir at low IL- 1 p. The therapy 
reconstitutes the existence of the pre-diabetic stable steady state 
(Fig. 7A, red circle). However, the original pathological steady state 
(Fig. 7A, black square) is neither deleted nor does it lose stability. The 
disease state remains stable irrespective of the glucose levels reached 
by glucose therapy. Therefore, the glucose-lowering therapy does not 
induce a transition from overt T2D to pre-diabetes and die functional 
P-ceU mass is further decreasing (Fig. 7D, black squares). 

The pre-diabetic state cannot be restored by IL-lRa 
therapy of overt T2D 

IL-lRa is a naturally occurring inhibitor of IL-ip, and it is 
upregnlated in obesity [28] and reduced in poorly controlled T2D 
patients [27]. By injecting exogenous IL-lRa, IL-lRa therapy aims at 
inhibiting IL- 1 P signalling. This is simulated in silico by fixing the level 
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Figure 6. A: The bifurcation diagram of the IL-ip-IL-1Ra subsystem as glucose (6) varies. This subsystem exhibits three stable (full lines) 
and two unstable (dotted blue lines) steady states. Above a critical glucose level (5.84 mM) the system loses stability of the pre-diabetic state (green 
line) and progresses to a high IL-1 (3 state (red line), vj/hich is stable across the whole physiological glucose range. A third stable state (black line) with 
even higher lL-1 (3 exists, which, in the hyperglycaemia range, may be associated with advanced T2D. B: A transition from five to three steady states is 
found that corresponds to the transition from pre-diabetes to overt T2D. Steady states appear as crossing points of the nullclines. The nullclines of IL- 
1 p (red) and IL-1 Ra (black) cross in five points at normal glucose level (4.5 mM). When glucose increases, the IL-1 p nullcline raises its minimum (inset, 
transition from red to blue to green curve) such that the low IL-1 p and low IL-1 Ra steady state vanishes. The system is forced to switch to a high IL-1 p 
state which may be associated with overt T2D. The corresponding glucose level is 4.5, 5.84 and 6.84 mM for the red, blue and green IL-1 p nullcline in 
the inset, respectively. 
doi:10.1371/journal.pcbi.1003798.g006 



of IL-lRa at a higher concentration (Fig. 7B, red cross). The simulated 
therapy slows down the loss of P-cell mass (Fig. 7D, red crosses) 
because increased IL-lRa drives IL-1|3 into a less pathological regime. 
The higher IL-lRa is set, the more IL-ip is reduced. At a threshold 
level of IL-lRa, which is determined by the peak of the IL-1 P nullcline 
(Fig. 6B, red Hne), the response of IL- 1 P becomes qualitatively different 
and it is reset below its physiological normal level (Fig. 7B, green 
triangle). Though tlie speed of P-cell loss is significantly reduced, it is 
not stopped (Fig. 7D, green triangles). Releasing IL-lRa fixation leads 



to a restoration of the T2D state associated with accelerated loss of P- 
cell mass (Fig. 7B, black square), because the stable pre-diabetic steady 
state associated with low IL- 1 P is not restored by this therapy (Fig. 7B), 

The pre-diabetic state can be restored by combined 
glucose-lowering and IL-lRa therapy 

A combined IL-lRa and glucose-lowering therapy (Fig. 7C) is 
required to induce functional P-cell mass reconstitution. Strong 
IL-lRa therapy induces a transition from the overt T2D steady 
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Figure 7. T2D therapy strategies and tlieir effect on p-cell mass. A: Glucose-lowering therapy. The pre-diabetic steady state is restored at the 
IL-ip (red line) and IL-IRa (black line) nullcline intersection (red circle). But the system stays in the stable T2D steady state (black square). B: 
Permanent IL-1 Ra therapy. The IL-1 p (red line) and IL-1 Ra (black dotted line) nullclines do not intersect at low IL-1 p and the system starts in the stable 
overt T2D steady state (black square). Clamping IL-IRa to levels below the peak of the IL-ip nullcline (lower blue line) reduces IL-ip (red cross). 
Clamping IL-1 Ra to levels above the peak (upper blue line) makes IL-1 p jump to low levels (green triangle). C: Combined IL-1 Ra and glucose-lowering 
therapy. Above threshold clamp of IL-IRa (dashed dotted line) induces a substantial IL-ip reduction (red dot). Simultaneous control of glucose 
restores a stable steady state at low IL-ip (red and black nullclines intersect at the blue circle). A release of IL-IRa clamp induces a transition to this 
steady state. D: Effect of the different treatments on p-cell mass. Therapy as in A (black squares), B (red crosses for below threshold IL-IRa clamp; 
green triangles for above threshold IL-1 Ra clamp), C (red dots for combined IL-1 Ra and glucose-lowering therapy; blue circles after release of IL-1 Ra 
clamp). No treatment (black line). 
doi:1 0.1 371 /journal.pcbi.1 003798.g007 



state to a transient state at low IL-ip (Fig. 7C, from black square 
to red dot) associated with reduced speed of P-ceU loss (Fig. 7D, 
red dots). A concurrent glucose-lowering therapy restores a stable 
steady state at low IL-1 (3 (Fig. 7C, blue circle). Upon withdrawal of 
the IL-IRa frxation, the in silico patient makes a transition to the 
restored stable steady state (Fig. 7C, from red dot to blue circle) 
which is associated with |3-ceU growth (Fig. 7D, blue circles), thus, 
inducing a long-standing improvement of the disease state. 

The IL-1 Ra therapy nnust be released to achieve a disease 
state improvement 

Long-term combined IL-IRa and glucose-lowering therapy 
prohibits P-cell mass regrowth in silico (Fig. 7D, red dots). The 
model predicts that IL- 1 Ra therapy, simulated in the model as the 
fixation of IL-IRa, must be released in order to allow for 
functional P-ceU mass reconstitution. P-ceU mass regrowth is 
suppressed during IL-IRa therapy as it reduces IL-ip to levels 
below the pre-diabetic state. 

Discussion 

An IL-ip switch model of T2D onset after long-term 
compensation was presented and a strategy for an anti-inflamma- 
tory therapy of T2D was proposed. The model is characterised by 



a pre-diabetic state during which increasing insulin resistance 
drives hyperglycaemia. This ultimately leads to a changed steady 
state configuration and progression into the overt T2D disease 
state. The model parameters are derived from in vitro islet 
measurements, from steady state conditions, and from glucose 
history of T2D patients, such that the model parameters are fuUy 
determined. In the following we first discuss its relevance to the 
empirical theory of five-stage P-ceU dysfunction [3] , then describe 
the predictions of the model, its relevance and potential 
contributions to clinical trials [23,26], and its relation to previous 
modelling works [33,35]. Last but not least, we discuss its 
limitations. 

Interestingly, the model shows three steady-states of IL-1 p. This 
result is beyond our expectations when designing the model but is 
consistent with the five-stage diabetes development theory [3]. 
This theory is mainly based on clinical observations, where stage 2, 
"adaptation" (termed as compensation in this paper) corresponds 
to the model state of lowest IL-ip (green in Fig. 6A). Stage 3 
"transient unstable early de-compensation" corresponds to the 
sudden switch of IL-ip/IL-lRa in the model, where a 55 fold 
increase in IL- 1 P level is accompanied with only 2 fold more IL- 
IRa, and the transient adaptation time thereafter. Stage 4 "stable 
de-compensation" corresponds to the state with medium IL-ip 
level (red in Fig. 6A) and stage 5 "severe de-compensation" 
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Figure 8. The dose-effect relationship of IL-1Ra therapy for 
mild (red circle) and strong T2D (black dot) in silico are shown. 

Following the protocol of the clinical trial, glycaemia improvement after 
90 days of therapy versus control, is used as a measure. The IL-IRa 
threshold predicted by the model is reflected in the jump of the 
resulting glucose level for both T2D groups. The optimal ILI-Ra dose for 
strong T2D is substantially larger because of the lower endogenous 
level of IL-IRa. The best achieved glucose improvement in response to 
the therapy is three-fold higher in the case of strong than in mild T2D. 
The starting glucose level is 11.96 and 18.02 mM for mild and strong 
T2D, respectively. The final glucose level of the both control groups is 
12.21 and 18.82 mlVl, respectively. 
doi:1 0.1 371 /journal.pcbi.1 003798.g008 

corresponds to the state with highest IL-ip level (black in Fig. 6A). 
Stage 1 "normal" is not represented in the model because glucose 
level is already higher even 14 years before T2D diagnosis (see 
Fig. 4A and [5]). Thus, the IL-1|3 switch model qualitatively 
reflects all stages of T2D development, starting from the 
compensation phase and progressing to the mild disease and 
ultimately to the strong disease state. 

Stage 3 is observed to be transient, persisting a few weeks in 
diflerent animal models that were subjected to partial pancreatec- 
tomies [41] or islets transplantation [42,43]. Interestingly, islets 
transplantation [43] and pancreatectomies [41] in rodents gave rise 
to glucose levels in two discrete groups: either normal glucose or 
severe diabetes, with almost no glucose levels in between, a pattern 
consistent with the idea of two distinct states in the presented model. 
Stage 3, transient and unstable, lies between the clinically non- 
pathologic and pathologic states, and thus, may be of great interest in 
terms of providing a time window during which therapies are more 
efficient in reversing the pathological processes right after diagnosis. 
Results from newly diagnosed diabetic patients who received 
intensive insulin therapy supports this notion [31]. However, in the 
current model, the switch oflL-lfi/IL-lRa, which is relevant to stage 
3, is instantaneous (see Fig. 6 and Fig. S2) because the model does not 
describe the dynamics of the IL- 1 fi/IL- 1 Ra subsystem during the 
transition. The equations of the IL-1 fi/IL- IRa are mainly based on 
assumed steady state relations, which were derived from in vitro data. 
In addition, the model does not describe the dynamics of insulin 
resistance, the improvement of which has been assumed to be 
responsible of preventing/postponing the transition during stage 3 
[3]. Rather, insulin resistance is implemented as an input to the 
glucose/insulin subsystem. These limitations of the current modeling 
approach (see also later discussions about the meaning of G and / in 
the model) make information about the duration of stage 3 
unavailable from the model. 



The IL- 1 fi bi-stable switch model suggests four requirements for 
a successful anti-inflammatory treatment of T2D: 

1 . Glucose-lowering therapy is combined with IL- 1 Ra therapy; 

2. IL-IRa is kept at an above threshold level; 

3. IL-IRa fixation is released when IL-lfi levels are down; 

4. Glucose-lowering therapy is continued on long-term. 

Then, the functional (i-ceU mass is predicted to increase. But 
even if it reaches its healthy level, the in silico T2D patient would 
remain in the pre-diabetic state, since the release of glucose control 
would restart T2D progression due to the pathological state of 
insulin resistance. The power of this therapy is not to cure T2D 
but to keep the in silico patient in the pre-diabetic state. 

Dose and duration of anti-IL-lpS therapy are critical. On the 
one hand, an above threshold IL- 1 Ra is required for inhibition of 
IL-ip auto-stimulation and the reverse transition from P-cell 
failure to compensation. On the other hand, IL- 1 fi has to be kept 
above its healthy level for P-ceU reconstitution. This infers a U- 
shaped dose-effect relationship of anti-IL-ip therapy in silico, 
which has been confirmed in a phase 1 clinical trial using a 
monoclonal antibody of human IL-ip [23]. 

Large scale clinical trials with anti-IL-lfS agents revealed that 
glycaemia was less efliciently reduced in patient groups with short (6 
years) compared to groups with long ( 1 1 years) disease duration 
[26]. We have performed simulations of IL-IRa therapy on patients 
with short and long T2D duration, where short and long T2D are 
associated with mild and strong T2D steady states, respectively 
(Fig. 8). The strong T2D state is initiated by manually increasing IL- 
1 fi and decreasing IL-IRa to the levels shown in Fig. 6B at 1 1 years. 
Then, anti-IL-ip therapy (ignoring other adjuvant medication) is 
applied by an IL-1 Ra fixation for 90 days to both in silico patient 
groups. The IL-IRa dose-efiect relationship for the two groups, 
expressed by the glycaemia improvement versus control, is 
compared (Fig. 8). In agreement with the clinical-trials results, 
optimal IL-IRa dosage induced more efficient glucose improve- 
ment for strong than for mild T2D. This happens because in strong 
T2D the aggravation of blood glucose levels is faster in the control 
group. The optimal dose of IL-IRa is different in both groups in 
silico (Fig. 8). This calls for an individualised treatment according to 
the endogenous inflammatory state. These simulations do not aim at 
explaining the diverse outcomes from clinical trials because other 
adjuvant medications are neglected during the simulations. 
Furthermore, the capability of the current model is limited to the 
low to medium hyperglycaemia/IL-ip region (see later discussion). 
However, the model stiU provides a framework in which different 
outcomes can be discussed in respect to patients' internal 
inflammatory levels and strategies may be suggested for patient 
stratification with respect to their inflammatory levels at the 
beginning of a clinical trial. 

The model proposed by De Gaetano et al [35] is able to account 
for different arms of the Diabetes Prevention Program (DPP) 
results [44] , by modeling the effects of different drugs on insulin 
resistance development. Some differences and similarities are 
worth noting here, besides those already discussed in the 
introduction. In both models glucose toxicity introduces a positive 
feedback in the disease progression. Therefore, both models 
emphasize the importance of an early treatment of hyperglycae- 
mia. The unique property of the current model is that the 
inflammatory signal, once switched on by glucose, cannot be 
switched off by glucose-lowering therapies. It has to be treated 
separately. This point seems to be consistent with recentiy 
discovered anti-inflammatory effects of some long-standing anti- 
diabetic agents. In fact, insulin has long been known for its anti- 
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inflammatory eflict, which, early presumed to be due to its 
glucose-lowering ability, recendy was established to be a comple- 
ment to its metabolic role [for a recent review see 45]. 
Additionally, there is emerging evidence suggesting that metfor- 
min, the first-line diabetic drug, have anti-inflammatory efTects 
that are independent of its hypoglycaemia effect [46], even though 
its anti-inflammatory effect in T2D being not as robust as in other 
diseases, such as hypertension [47^9] . The impact of currendy 
employed T2D drugs onto the inflammatory state of the patients 
prohibits falsification or verification of the model prediction that 
the T2D disease state cannot be alleviated by controlling 
glycaemia alone. 

The presence of the IL-ip switch suggests that the system is 
operating in two distinct regimes, which defines a clear threshold 
for an anti-inflammatory therapy. A striking prediction of the 
model is that the anti-inflammatory therapy has to be released 
after a certain time. These inferences are derived from the IL-lfi/ 
IL-lRa hysteresis and are consistent with some long-standing 
findings [11,12], as well as recent progresses in the field, such as 
multistage development of the disease [3,5], long-term effect of 
temporary IL-lRa in overt-diabetes [20] and the efficiency of 
glucose control in pre-diabetes but not overt-diabetes [31,36]. 

The present IL-ip switch model is built on in vitro data and 
several assumptions, both of which cast limitations on its 
applicability. The in vitro data show regulation of IL-1(3 and IL- 
lRa by glucose and exogenously added IL-ip [14,17,27,32]. 
Other factors involved in in vivo IL-ip control, such as free fatty 
acids [50], infiltrated macrophages and islet amyloid polypeptides 
[51] are not considered. These factors are associated with serious 
hyperglycaemia and overt T2D. However, they do not or only 
weakly affect the results of the core IL-ip model, especially 
concerning the onset bifurcation at glucose levels below 6 mM. 

IL- 1 P induced P-cell proliferation and apoptosis were extrap- 
olated from data of cultured islets [32] . It is unclear to what degree 
the in vitro data can reflect relevant in vivo processes, but several 
weak points of this assumption could be delineated. There is 
evidence that P-cells can re-enter cell cycle in vivo and that P-ceU 
mass can adapt to metabolic needs [52], both supporting adaptive 
P-ceU proliferation in vivo. However, it is not the only possible 
source of functional P-cells. Dediflferentiation of other cell types 
with a common precursor [53], as well as reactivation of 
dysfunctional P-cells [8], contributes to P-ceU homeostasis. 
Furthermore, data from autopsy indicate that the proliferation 
rate of existing P-cells is very low. Consequentiy, increased P-ceU 
mass seen in the compensation phase may be mainly attributed to 
new islet formation from other cell types [54]. A more realistic 
description of P-ceU mass dynamics may further improve the 
model and extend its vaUcUty to later stages of the disease. 

By representing the inflammatory stimulus via F instead of IL-1 P, 
it is implicidy assumed that the amount of membrane IL- 1 P receptor 
per P-ceU is constant. In support of this, there- is c\ idence that the 
amount of IL-1 P receptor in the P-ceU membrane is much higher 
than in other cell types [55], suggesting that it is saturated. However, 
it is unlikely that the receptor is not down-regulated at high 
concentrations of IL-1 p. Consequentiy, the model may overestimate 
the inflammatory signal at high IL- 1 P levels. This limitation has to be 
kept in mind when interpreting model results at high IL- 1 p. 

Another important assumption concerns the IL-ip mRNA 
quantity used in the model. The IL- 1 P mRNA quantit)- is assumed 
to be proportional to the IL- 1 P protein level, in order to establish the 
quantitative relationships in IL-ip production. However, IL-ip 
production is characterized by the dissociation of transcription and 
translation. As a result, high levels of IL-ip mRNA are always 
associated with low levels of protein [29]. Alternatively, a sigmoidal 



relationship may be considered instead of the linear one used here. 
While it is difficult to estimate the quantitative error induced by the 
assumed linear relationship, it is straightforward to predict that the 
alternative sigmoidal relationship would increase the nonlinearit^- of 
IL- 1 P auto-stimulation, and of the interaction between glucose and 
IL-1 p. Therefore, the essential elements for the bi-stable switch 
remain unchanged (or even would be increased), such that the 
quaUtative results of the IL-ip switch model are robust. 

It is difficult to define the exact physiological meaning of G and 
/ in the model [35]. Since G and / are described by one- 
compartment dynamics in the model, only changes in long-term 
fasting levels are accessible to the model. Hence, we propose to 
interpret G as an abstract parameter associated with long-term 
fasting glucose on a time scale of days to weeks. Fasting levels 
represent more the internal property of the body's glucose control 
system than the history of how it was perturbed by meal or 
exercise [56]. However, the driving stimulus of IL-1 P production 
may depend on peaks of glycaemia rather than on fasting levels. It 
is possible that dysregulations of fasting glucose levels and 
postprandial peaks are independent, as evidenced by patients 
with impaired glucose tolerance but normal fasting glucose. 
Eventually, the modifications of the dynamics of insulin release 
(first-phase insulin) are more important for glycaemia control than 
the long-term fasting levels of glucose. We would like to emphasize 
that these detaUs of glucose and insuUn dynamics are beyond the 
resolution of the model presented here and have to be addressed 
with a model focusing on the short-term dynamics of these 
quantities. 

T2D is a very complex disease. The current model is simplified 
on purpose in order to generate structural insight. Other players of 
the glucose control system, for example glucagon, also play very 
important roles in glucose homeostasis and T2D development. 
Augmented glucagon secretion, together with impaired P-ceU 
function, happens in the very early phase of pre-diabetes [57,58]. 
The underlying islet defects were identified as reduced maximal 
insulin response and reduced glucose-sensitivity of P-/a-ceU [58]. 
According to the physiological integral rein control theory [56,59], 
glucose-sensitivity of P-/a-ceU are important factors in setting 
steady state glucose level. In fact, recent studies support an 
independent role of glucose-sensitivity in hyperglycaemia devel- 
opment [60,61]. Future improvement of the model should include 
the glucose-sensitivity of P-/a-cell to account for the rising of 
fasting glucose in the compensation phase. In particular, this might 
improve the fitting to the data in Fig. 4A during the compensation 
phase. 

In addition, post-load glucose of incident diabetes cases showed 
a rapid increase 5 years before diagnosis [5], which is followed by 
a rapid decrease in insulin sensitivity. This important change is 
not reflected in the history of fasting glucose whUe it is 
represented in the model output (Fig. 4A). On the one hand, 
the absence of a corresponding change in fasting glucose at ~5 
years before diagnosis suggests that fasting glucose level is 
dominated by factors other than insulin resistance. Again, this 
might point to a role of glucagon. On the other hand, the 
coincident post-load glucose increase and insulin sensitivity 
decrease at ~5 years before diagnosis might suggest a cause- 
effect relationship between P-ceU function and insulin resistance 
[62] . The complexity of the defects during the pre-diabetes phase 
highlights the necessity of the mathematical modeUng approach 
in the field. 

Although results from both in vitro experiments and animal 

models were promising, the actual effect of anti-IL- 1 P agents in 
clinical trials of T2D was a matter of debate because the observed 
effect was very modest so far. Besides the possible role of 
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polymorphisms of the ILIRN gene which encodes IL-lRa [20], 
more research is needed to clarify for example, drug interaction 
with other IL-1 family members, such as soluble IL-1 receptors 
and soluble IL-1 receptor accessory proteins, in the pancreatic 
local environment as well as in blood and liver. 

On the structural level, the present analysis proposes that the 
transition from pre-diabetes to overt T2D is associated with 
changed stability properties of the endocrine-immune system. 
The transition is initiated at a threshold value of glucose and 
associated with qualitatively different inflammatory states. The 
notion of an IL- 1 fi hysteresis is supported by data from different 
studies, which showed the effectiveness of glucose control in pre- 
diabetes [36] or newly diagnosed patients [31] and the 
ineffectiveness of glucose control in long-standing T2D patients 
[11-13]. It is also supported by the clinical trial based on a 
combined IL-lRa and glucose-lowering therapy: the improve- 
ment of the disease state remained for a year after release of the 
IL-lRa therapy [20]. The existence of the IL-ip hysteresis 
requires a combined glucose-lowering and anti-inflammatory 
therapy for overt T2D patients, where the anti-inflammatory 
therapy should be kept for only a limited time. 
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